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^ Abstract 

^ This paper, together with a subsequent paper, questions the so-called 'LSND anomaly' : 

^ a 3.8 (T excess of Uc interactions over standard backgrounds, observed by the LSND Col- 

^ laboration in a beam dump experiment with 800 MeV protons. That excess has been inter- 

w preted as evidence for the — ^ z^c oscillation in the Am? range from 0.2 eV^ to 2 eV^. 

^ Such a Am^ range is incompatible with the widely accepted model of oscillations between 

^ three light neutrino species and would require the existence of at least one light 'sterile' neu- 

I— I trino. In this paper, new data on pion production by protons on nuclei are presented, and 

^ four decades old data on pion production by neutrons on nuclei are recalled, that together 

I increase significantly the estimates of standard backgrounds in the LSND experiment, and 

^ decrease the significance of the 'LSND anomaly' from 3.8 a to 2.9 a. In a subsequent 

4^ paper, in addition the LSND Collaboration's data analysis will be questioned, rendering a 
further reduction of the significance of the 'LSND anomaly' . 
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1 Introduction 

In the years 1993-1998, the LSND experiment was carried out at the Los Alamos National 
Laboratory. Its scientific goal was a search for z/^ — )■ Uc oscillations in the 'appearance' mode. 
The neutrino fluxes were produced by dumping 800 MeV protons into a 'beam stop'. While z/^, 
and Uc fluxes were abundant, the Ue flux was vanishingly small. 

The LSND Collaboration claimed an excess of interactions over the expectation from 
standard backgrounds [1]. This excess was interpreted as evidence for the z/^ — )■ z>c oscillation 
with Am^ in the range from 0.2 eV^ to 2 eV^ and came to be known as 'LSND anomaly'. In 
stark conflict with the widely accepted model of oscillations of three light neutrino flavours, the 
excess would require the existence of at least one light 'sterile' neutrino that does not couple to 
the Z boson. 

Since the 'LSND anomaly' calls the Standard Model of particle physics in a non-trivial way 
into question, the MiniBooNE experiment at the Fermi National Accelerator Laboratory set out 
to check this result. In the neutrino mode, an oscillation — )• with parameters compatible 
with the LSND claim was not seen [|2|. However, first results from running in antineutrino mode 
and searching for z/^ — > Uc led the MiniBooNE Collaboration to conclude that their result does 
not rule out the 'LSND anomaly' [3] that had indeed been observed in the antineutrino mode. 

In this situation it appears worthwhile to undertake a critical review of the original results 
of the LSND experiment that gave rise to the 'LSND anomaly'. This is the subject both of this 
paper and of a subsequent paper Q. 

Our review of the physics and technical papers published by the LSND Collaboration led to 
agreement with many of their approaches and results. However, in two areas we disagree with 
LSND. The first area — which is discussed in this paper — concerns the underestimation of the 
z/o flux from standard sources. The second area — which will be discussed in Ref. p| — concerns 
the underestimation of systematic errors in the isolation of the signal of ~120 z/e + p — e+ + n 
reactions with a correlated 7 from neutron capture out of ~2100 candidate events. 

Our review of the LSND claims was carried out without insider information from members 
of the LSND Collaboration^^ We came to consider this advantageous because it led us to 
reproduce and assess from an entirely independent point of view the steps that were taken in the 
LSND analysis. 

In view of the importance of the LSND result for neutrino physics, the Leitmotiv of our 
work was to answer the question: did LSND provide "Evidence for neutrino oscillations from 
the observation of z^c appearance in a z/^ beam" [JJ? 

Our short answer is no. Our long answer is laid down in this and a subsequent paper. 

2 Neutrino fluxes in the LSND experiment 

In order to search for the oscillation z/^j — )• z/c in the 'appearance' mode, a negligible flux of z/g 
is needed amidst a large flux of z/^. The LSND experiment was designed to realize this situation 
as closely as possible by dumping a large number of protons of 800 MeV kinetic energy into a 
totally absorbing 'beam stop' made essentially of water, copper and iron. 

The interaction of the beam protons with (primarily) O, Cu and Fe nuclei led to a hadronic 
cascade, schematically depicted in Fig. [T| The relatively low beam proton energy of 800 MeV 
assured that the hadrons in the cascade were only protons, neutrons, 7r+ and tt^, and that 

^'Excepting the provision by M. Sung of a FORTRAN program with the LSND parametrization of pion produc- 
tion cross-sections, and of the LSND Technical Note 134 (Neutrino Fluxes for LSND, by Myungkee Sung, Rex 
Tayloe and David Works, May 2000). 
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Fig. 1: Schematic display of tlie 1st and 2nd generation hadrons in the cascade after the interaction of an 
800 MeV beam proton with a nucleus in the LSND 'beam stop'. 
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Fig. 2: Schematic display of the decay processes of positive (upper chain) and negative (lower chain) 
pions that are created in the hadronic cascade after the interaction of an 800 MeV beam proton with a 
nucleus in the LSND 'beam stop'. 



the number of 7r+ was an order of magnitude larger than the number of vr^ . The totally ab- 
sorbing target assured that the flight path of pions — the only short-lived particles among those 
hadrons — was short so that there was a small probability to decay in flight or to interact, but a 
large probability to come to rest after their kinetic energy was dissipated by ionization. 

The fate of charged pions in the hadron cascade after the interaction of the beam proton are 
schematically shown in Fig. [2j Since there are decisive differences between what happens to 
positive and negative pions and their decay products, we discuss them in turn. 

Positive pions (the upper chain in Fig. [2]) come to rest and diffuse around until their decay 
Ti~^ — i- This creates an isotropic flux of monoenergetic z/^ of 30 MeV and monoenergetic 

fi^. Because of the relatively long lifetime of muons, essentially all /i+ will come to rest after 
their kinetic energy is dissipated by ionization. They also will diffuse around until their decay 
yLi+ — >■ e'^u^i'^. This creates isotropic fluxes of and with well known respective energy 
spectra. Thus 7r+ lead to large fluxes of u^, and u^. 

By contrast, negative pions (the lower chain in Fig. |2]) that come to rest, are caught in low- 
lying atomic orbits and disappear by strong interaction with the respective atomic nucleus, with 
no neutrino emitted. 

Additional, albeit much less abundant, neutrino fluxes stem from pions and muons that de- 
cay in flight before they interact or come to rest. For positive pions and muons, that means only 
a small addition to the fluxes of v^, Vc and from decays at rest. Of much larger importance, 
however, are the analogous fluxes of v^, Uc and from the decays of negative pions and muons 
where the negative pion decayed in flight. 

Neutrinos from negative muon decays can either come from decays in flight or from decays 
at rest. In the latter case, muon decay competes with the weak interaction of muons that are 
caught in low-lying atomic orbits with the respective nucleus, /i~ -i- (A, Z) ^ + {A, Z — 1), 
that leads to an isotropic z/^j flux with an energy around 100 MeV which adds to the — already 
abundant — flux of and therefore is of secondary interest. Of primary interest, however, is 
the isotropic flux of z/c that arises from the decay of those /i^ in orbit that escape the capture 
reaction with the nucleus. This small flux constitutes a conventional background to the search 
for the — )■ oscillation and therefore must be well understood^\ 

The conventional flux crucially depends on 
^''In comparison, Vc from the decays of negative muons in flight are, although taken into account in our analysis. 
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1 . the total number of tt in the hadron cascade; 

2. the momentum spectra of vr" (the decay probability per unit length is inversely propor- 
tional to the momentum); 

3. the angular spectra of tt" (different drift spaces and materials are sampled at different 
angles); and 

4. the geometrical layout and the material composition of the 'beam stop' (in free drift 
spaces there is decay but virtually no interaction and no dissipation of energy by ion- 
ization; the competition between capture and decay of negative muons depends strongly 
on the atomic number Z). 

As for hadron production properties, the challenge is quite daring: for protons, neutrons, and 
charged pions, the interaction cross-sections and the double-differential inclusive cross-sections 
of secondary hadron production are needed to calculate neutrino fluxes, in the kinetic energy 
range between 800 and zero MeV, for several 'beam stop' materials, primarily water, copper 
and iron. 

At the time of the LSND experiment, pertinent data were scarce. 

At CERN, it was realized that the HARP experiment could make a significant contribution 
to items 1 to 3 of the above list, by measuring pion production in the interactions of 800 MeV 
(=1500 MeV/c) protons with relevant target nuclei^\ Being not part of the experiment's ap- 
proved programme, only a few days of data taking could be set aside for the study of this issue. 

In Section[3]we discuss briefly those features of the HARP experiment that were used to ana- 
lyze the pertinent data. In Section|4[ we report on the cross-sections that the HARP-CDP group 
extracted from these data, and that are used to cross-check the calculations of the conventional 
background published by LSND [T|. 

3 The HARP experiment 

The HARP experiment was designed to carry out a programme of systematic and precise (i.e., at 
the few per cent level) measurements of hadron production by protons and pions with momenta 
from 1 .5 to 15 GeV/c, on a variety of target nuclei. It took data at the CERN Proton Synchrotron 
in 2001 and 2002. Cross-sections of large-angle pion, proton and deuteron production on a 
number of nuclei have been published (see Refs. [SmTll). 

For the work described here, a thin and a thick water target was exposed to a +1.5 GeV/c 
beam, with length of 6 cm and 60 cm, respectively. The target radius was in both cases 18 mm. 
The water was contained in plexiglass tubes with 2 mm wall thickness. The thin water target 
was used for the measurement of the double-differential inclusive pion production cross-section 
by protons. The finite thickness of the thin water target leads to a small attenuation of the 
number of incident beam particles. The attenuation factor taken into account in the cross-section 
normalization is /att = 0.964. 

The thick water target was used for a relative measurement: a high- statistics determination 
of the 7r~/7r"'" ratio of secondary pions. 

In addition to the water targets, data were recorded with the same thin (5% Aint) tantalum, 
copper and lead targets as used for the work described in Refs. [|7]-[9). 

Protons and pions with 1.5 GeV/c momentum were delivered by the T9 beam line in the 
East Hall of CERN's Proton Synchrotron with a momentum bite Ap/p ~ 1%. Only the pos- 
itive beam polarity was used. The beam instrumentation, the definition of the beam particle 

nearly negligible. 

"^^We thank G.B. Mills from the Los Alamos National Laboratory for his suggestion to take such data. 
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trajectory, the cuts to select 'good' beam particles, and the muon and electron contaminations 
of the particle beams, are as described, e.g., in Ref. [5J. 

Our calibration work on the HARP TPC and RPCs is described in detail in Refs. []T2| and 



[ 13 1, and in references cited therein. The momentum resolution cr(l/j9T) of the TPC is typically 
0.2 (GeV/c)"^ and worsens towards small relative particle velocity (3 and small polar angle 9. 
The absolute momentum scale is determined to be correct to better than 2%, both for positively 
and negatively charged particles. 

The polar angle 9 is measured in the TPC with a resolution of ~9 mrad, for a representative 
angle of ^ = 60°. In addition, a multiple scattering error must be considered that is for a proton 
with pt = 500 MeV/c in the TPC gas ~4 mrad at 6* = 20°, and ~12 mrad at 6* = 90°; for 
a pion with the same characteristics, the multiple scattering errors are ~3 mrad and ~6 mrad, 
respectively (the numbers refer to the thin water target). The polar-angle scale is correct to 
better than 2 mrad. The TPC measures dE/dx with a resolution of 16% for a track length of 
300 mm. 

The efficiency of the RPCs that surround the TPC is better than 98%. The intrinsic time 
resolution of the RPCs is 127 ps and the system time-of-flight resolution (that includes the jitter 
of the arrival time of the beam particle at the target) is 175 ps. 

To separate measured particles into species, we assign on the basis of dE/dx and (3 to each 
particle a probability of being a proton, a pion (muon), or an electron, respectively. The proba- 
bilities add up to unity, so that the number of particles is conserved. These probabilities are used 
for weighting when entering tracks into plots or tables. Because of the strong preponderance of 
protons over 7r+, great care had to be taken to separate cleanly these two particles. To illustrate 
the difficulty of this undertaking. Fig. [3] shows the statistics of secondary protons, 7r+, vr^, e+ 
and e^, with polar angles in the range 20° < 9 < 60°, for beam protons with 800 MeV kinetic 
energy impinging on the thin water target. 

4 Pion production by 1.5 GeV/c protons 
4.1 Pion cross-sections on various target nuclei 

In Tables [I]-[8} we give the double-differential inclusive cross-sections d^a/dpdfi of 7r+ and vr^ 
production in thin targets of water, copper, tantalum, and lead, for incoming beam protons and 
beam 7r+, including statistical and systematic errors. The cross-sections are given in bins of pq^, 
for two ranges of polar angle: 20° < 9 < 60° and 60° < 9 < 125°. The cross-sections for the 
polar-angle range 20° < 9 < 60° are graphically shown in Fig.|4j 

Because of the limited time available for data taking, the statistical precision of the re- 
ported cross-sections is rather poor. Systematic errors stem primarily from uncertainties in the 
momentum scale and the absorption of secondary pions in materials between the vertex and the 
sensitive TPC volume. Smaller systematic errors stem from the momentum resolution, the time- 
of-flight and dE/dx resolutions, the subtraction of secondary protons and electrons, the sub- 
traction of beam muons (incoming 7r+ beam only), and overall normalizations. All systematic 
errors are propagated into the spectra of secondary pions and then added in quadrature. They 
add up to a systematic uncertainty of our inclusive cross-sections of ~7% for 20° < 9 < 60°, 
and of ~1 1% for 60° < ^ < 125°, respectively. 

The question that is most interesting in the context of the 'LSND anomaly' is a comparison 
of our measured HARP-CDP cross-sections with those that were used by LSND in the calcu- 
lation of the neutrino fluxes in their experiment. There was a considerable effort to take into 
account all relevant data that existed at that time, and to merge them into what was considered 



by LSND a reliable parametrization of pion production by protons on nuclei [ 14- 17 1. Their 
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Fig. 3: Statistics of secondary protons, 7r+, tt , e+ and e for beam protons with 800 MeV kinetic 
energy impinging on water. 
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Fig. 4: Double-differential inclusive cross-section d^a/dpdO. [mb/(GeV/c sr)] of the production of vr^ 
by incoming 1 .5 GeV/c 7r+ and protons impinging on thin water (upper left panel), copper (upper right 
panel), tantalum (lower left panel) and lead (lower right panel) targets. 
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Table 1: Double-differential inclusive cross-section d^cj/dpdri [mb/(GeV/c sr)] of the production of 
pions in p -I- H2O — vr^ -1- X interactions with 1.5 GeV/c beam momentum; the first error is statistical, 
the second systematic; pt in GeV/c, polar angle 9 in degrees. 





20 < e < 60 


Pt 


7r+ 




0.10-0.15 
0.15 - 0.20 
0.20-0.25 
0.25 - 0.35 


33.6 ± 1.9 ±2.1 

53.7 ±2.3 ±3.3 
47.2 ± 2.1 ± 2.9 
15.4 ±0.9 ±0.9 


7.0 ± 0.9 ± 0.4 
6.7 ± 0.8 ± 0.4 
3.5 ± 0.6 ± 0.2 

1.1 ±0.2 ±0.1 




60 < < 125 


PT 




7r 


0.10-0.15 
0.15 - 0.20 
0.20-0.25 
0.25 - 0.35 


23.9 ± 1.3 ± 1.5 

32.4 ± 1.4 ±2.0 

21.5 ± 1.2 ± 1.3 
7.5 ± 0.5 ± 0.5 


5.6 ± 0.7 ± 0.3 
5.4 ± 0.6 ± 0.3 
3.2 ± 0.4 ± 0.2 
0.9 ±0.2 ±0.1 



Table 2: Double-differential inclusive cross-section d^cj/dpd$7 [mb/(GeV/c sr)] of the production of 
pions in 7r+ -1- H2O tt^ +X interactions with 1.5 GeV/c beam momentum; the first error is statistical, 
the second systematic; pT in GeV/c, polar angle 9 in degrees. 





20 < e < 60 


PT 


n+ 


n 


0.10 - 0.15 
0.15 - 0.20 
0.20 - 0.25 
0.25 - 0.35 
0.35 - 0.45 


51.8 ± 1.7 ± 3.3 
70.1 ± 1.8 ±4.5 

71.8 ± 1.8 ±4.6 
59.4 ± 1.2 ±3.8 

50.9 ± 1.2 ±3.3 


29.2 ± 1.3 ± 1.9 
26.8 ± 1.1 ± 1.7 

26.7 ± 1.1 ± 1.7 
22.5 ± 0.7 ± 1.5 

13.8 ±0.6 ±0.9 




60 < e < 125 


PT 


n+ 




0.10 - 0.15 
0.15 - 0.20 
0.20 - 0.25 
0.25 - 0.35 
0.35 - 0.45 


43.4 ± 1.3 ±2.8 
50.0 ± 1.3 ±3.2 

39.0 ± 1.1 ±2.5 

24.1 ±0.6 ± 1.6 
18.9 ± 0.6 ± 1.2 


23.1 ± 1.0 ± 1.5 
22.1 ±0.8± 1.4 
16.3 ± 0.7 ± 1.1 
10.0 ± 0.4 ± 0.6 
5.1 ± 0.3 ± 0.3 



parametrization is compared in Fig.|5]with our measurements. Given the notorious uncertainty 
of hadronic cross- sections especially in this low-energy domain, one might argue that the agree- 
ment is astonishingly good. In more detail, there are worrisome discrepancies visible, especially 
for heavier target nuclei, where the all-important tt production exceeds the one assumed by 
LSND. 

4.2 TT /tt^ ratio from proton interactions in water 

Because of the importance of the tt production rate in the interactions of protons with 800 MeV 
kinetic energy in the LSND 'beam stop' for the prediction of the conventional flux, in 
the HARP experiment a high-statistics measurement of the tx^ jix^ ratio in the interactions of 
1 .5 GeV/c protons with a 60 cm thick water target was undertaken. A water target was cho- 
sen since for most of the LSND data taking a thick water target contributed significantly to 
interactions of beam protons within the 'beam stop', as will be detailed in Section[6j 

Table|9]lists the ratios of inclusive double-differential cross-sections d^cr/dpdi7 of tt^ to 7r+ 
production by 1.5 GeV/c protons in the thick water target. Because of the thickness of the target, 
reinteractions of secondaries occur and contribute to the observed pion yields. These ratios are 
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Fig. 5: Comparison of measured HARP-CDP 7r+ (black points) and vr^ (open points) production cross- 
sections on water, copper, tantalum and lead, with the ones predicted by the LSND parametrization (full 
lines for 7r+, broken ones for 7r~), at the polar angles of 40° and 90°. 
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Table 3: Double-differential inclusive cross-section d^d/dpdri [mb/(GeV/c sr)] of the production of 
pions in p -I- Cu — vr^ -I- X interactions with 1.5 GeV/c beam momentum; the first error is statistical, the 
second systematic; in GeV/c, polar angle 9 in degrees. 





20 < e < 60 


Pt 


71+ 


n 


0.10 - 0.15 
0.15 - 0.20 
0.20 - 0.25 
0.25 - 0.35 


43.8 ± 5.9 ± 3.2 
53.2 ±5.8 ±3.9 
42.1 ±5.2 ±3.1 
19.8 ± 2.5 ± 1.4 


22.8 ± 4.4 ± 1.7 
20.2 ± 3.6 ± 1.5 
8.2 ± 2.3 ± 0.6 
1.7 ±0.7 ±0.1 




60 < 6» < 125 


PT 


n+ 




0.15 - 0.20 
0.20 ~ 0.25 
0.25 - 0.35 


36.2 ± 4.0 ± 4.0 
22.5 ± 2.9 ± 2.5 
8.6 ± 1.2 ± 1.0 


10.5 ± 2.1 ± 1.2 
4.7 ± 1.3 ±0.5 
2.3 ± 0.6 ± 0.3 



Table 4: Double-differential inclusive cross-section d'^a/dpdO, [mb/(GeV/c sr)] of the production of 
pions in 7r+ -I- Cu — )• vr^ -I- X interactions with 1.5 GeV/c beam momentum; the first error is statistical, 
the second systematic; in GeV/c, polar angle 9 in degrees. 





20 < < 60 


PT 


7r+ 


7r 


0.10 - 0.15 
0.15 - 0.20 
0.20 - 0.25 
0.25 - 0.35 
0.35 - 0.45 


96.2 ± 6.0 ± 7.3 
1 10.0 ± 5.9 ± 8.3 

96.3 ± 5.3 ± 7.3 
86.7 ± 3.5 ± 6.5 
78.7 ±3.6 ±5.9 


58.9 ± 4.6 ± 4.4 
54.0 ±4.0 ±4.1 

49.0 ± 3.8 ± 3.7 

37.1 ±2.3 ±2.8 
27.7 ± 2.1 ± 2.1 




60 < < 125 


PT 


71 + 


7r 


0.15 - 0.20 
0.20 - 0.25 
0.25 - 0.35 
0.35 - 0.45 


92.5 ± 4.3 ± 10.5 
74.5 ± 3.7 ± 8.5 
41.8 ± 1.9 ±4.7 
29.8 ± 1.7 ± 3.4 


54.2 ± 3.3 ± 6.2 

40.0 ± 2.7 ± 4.5 
21.8 ± 1.4 ±2.5 

13.1 ± 1.1 ± 1.5 



deemed to be useful for the simulation of the hadron cascade in the LSND 'beam stop' even if 
the measured tt^ to 7r+ production ratio does not exactly reproduce the situation in the LSND 
'beam stop', for differences in the water target geometries. 

Figure |6] shows graphically the ratios listed in Table [9]for the polar-angle range 35° < 9 < 
50°. 
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Fig. 6: Ratios of inclusive double-differential cross-sections d^cr/dpciJl of tt to tt"^ production by 
1.5 GeV/c protons in the thick water target, in the polar-angle range 35° < 9 < 50°. 
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Table 5: Double-differential inclusive cross-section d'^a/dpdfl [mb/(GeV/c sr)] of the production of 
pions in p + Ta — > TT^ + X interactions with 1.5 GeV/c beam momentum; the first error is statistical, the 
second systematic; pt in GeV/c, polar angle 9 in degrees. 





20 < 6* < 60 


Pt 


7r+ 


TT 


0.10 - 0.15 
0.15 - 0.20 
0.20 - 0.25 
0.25 - 0.35 


60.6 ± 10.4 ± 4.4 
70.6 ± 10.5 ± 5.1 
63.1 ±9.8 ±4.6 
24.9 ± 4.3 ± 1.8 


26.3 ± 6.6 ± 1.9 
33.9 ± 7.3 ± 2.5 
16.7 ± 5.1 ± 1.2 
5.2 ± 2.0 ± 0.4 




60 < < 125 


Pt 


7r+ 


TT 


0.15 - 0.20 
0.20 - 0.25 
0.25 - 0.35 


46.0 ± 7.4 ± 5.1 
34.6 ± 5.9 ± 3.9 
12.2 ± 2.3 ± 1.4 


23.3 ± 5.4 ± 2.6 
11.3 ±3.4 ±1.3 
3.3 ± 1.2 ± 0.4 



Table 6: Double-differential inclusive cross-section d'^a/dpdO, [mb/(GeV/c sr)] of the production of 
pions in 7r+ -I- Ta — )■ vr^ -I- X interactions with 1.5 GeV/c beam momentum; the first error is statistical, 
the second systematic; pt in GeV/c, polar angle 9 in degrees. 





20 < 6» < 60 


Pt 


77+ 


77 


0.10 - 0.15 
0.15 - 0.20 
0.20-0.25 
0.25 - 0.35 
0.35 - 0.45 


138.3 ± 11.9 ± 10.4 

137.2 ± 10.2 ± 10.4 

132.3 ± 9.6 ± 10.0 
114.8 ±6.4 ±8.7 
117.7 ±6.8 ±8.9 


98.1 ± 9.5 ± 7.4 
93.1 ±8.3 ±7.0 
68.3 ± 6.9 ± 5.2 
65.8 ± 4.9 ± 5.0 
48.3 ± 4.4 ± 3.6 




60 < 9 < 125 


Pt 


77+ 


77 


0.15 - 0.20 
0.20 - 0.25 
0.25 - 0.35 
0.35 - 0.45 


142.0 ± 8.9 ± 16.1 
123.4 ± 7.8 ± 14.0 
67.6 ± 4.0 ± 7.7 
45.9 ± 3.4 ± 5.2 


93.8 ± 7.6 ± 10.7 
56.4 ± 5.1 ± 6.4 

35.4 ± 2.9 ± 4.0 

22.5 ± 2.3 ± 2.6 



Table 7: Double-differential inclusive cross-section d^cr/dpdJl [mb/(GeV/c sr)] of the production of 
pions in p + Pb — )• 7r^ + X interactions with 1.5 GeV/c beam momentum; the first error is statistical, the 
second systematic; pt in GeV/c, polar angle 9 in degrees. 





20 < < 60 


PT 


77 + 


77 


0.10 - 0.15 
0.15 - 0.20 
0.20 - 0.25 
0.25 - 0.35 


49.1 ± 9.8 ± 3.6 

54.2 ±9.6 ±3.9 
56.0 ± 10.2 ± 4.1 
23.2 ± 4.6 ± 1.7 


44.5 ± 9.8 ± 3.2 
35.3 ± 7.9 ± 2.6 
12.9 ± 4.9 ± 0.9 
6.6 ± 2.3 ± 0.5 




60 < < 125 


Pt 


77 + 


77 


0.15 - 0.20 
0.20 - 0.25 
0.25 - 0.35 


41.4 ±7.2 ±4.6 

43.4 ± 6.8 ± 4.9 

13.5 ± 2.4 ± 1.5 


30.5 ± 6.1 ± 3.4 
9.0 ± 3.0 ± 1.0 
4.2 ± 1.4 ± 0.5 
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Table 8: Double-differential inclusive cross-section d'^a/dpdfl [mb/(GeV/c sr)] of the production of 
pions in 7r+ + Pb ^ 7r='= + X interactions with 1.5 GeV/c beam momentum; the first error is statistical, 
the second systematic; pr in GeV/c, polar angle 6 in degrees. 





20 < < 60 


Pt 


7r+ 




0.10 - 0.15 
0.15 - 0.20 
0.20 - 0.25 
0.25 - 0.35 
0.35 - 0.45 


134.8 ±11.9 ±10.2 
136.2 ± 10.2 ± 10.3 
152.4 ± 11.0 ± 11.5 
125.4 ± 6.9 ± 9.5 
131.2 ±7.7 ±9.9 


100.4 ± 10.1 ± 7.6 
81.9 ±8.1 ±6.2 
72.9 ± 7.5 ± 5.5 
68.7 ±5.1 ±5.2 
58.2 ± 5.0 ± 4.4 




60 <9 < 125 


Pt 


7r+ 


TT 


0.15 - 0.20 
0.20 - 0.25 
0.25 - 0.35 
0.35 - 0.45 


149.0 ± 9.5 ± 16.9 
104.9 ±7.3 ±11.9 
69.6 ±4.1 ±7.9 
53.9 ± 3.7 ± 6.1 


96.9 ±7.5 ±11.0 
81.6 ± 6.5 ± 9.3 
40.0 ± 3.2 ± 4.5 
25.3 ± 2.5 ± 2.9 



Table 9: Ratio of inclusive double-differential cross-sections d'^a/dpdU of tt~ to 7r+ production by 
1.5 GeV/c protons in a 60 cm thick water target; the first error is statistical, the second systematic; pT in 
GeV/c, polar angle 6 in degrees, fiducial z range in target in cm. 



Pt 


20 <e < 35 and 3 < z < 56 


35 < e < 50 and 3 < z < 59 


0.100-0.125 
0.125-0.150 
0.150-0.175 
0.175 - 0.200 
0.200 - 0.225 
0.225 - 0.250 
0.250 - 0.275 
0.275 - 0.300 


0.203 ± 0.018 ± 0.008 
0.112 ±0.010 ±0.007 
0.119 ±0.010 ±0.005 
0.065 ± 0.007 ± 0.004 
0.041 ± 0.006 ± 0.002 
0.064 ±0.011 ±0.003 
0.074 ± 0.015 ± 0.010 
0.076 ± 0.021 ± 0.013 


0.270 ± 0.028 ± 0.023 
0.210 ± 0.018 ± 0.008 
0.169 ±0.014 ±0.010 
0.126 ±0.011 ±0.005 
0.101 ±0.010 ±0.007 
0.072 ± 0.008 ± 0.003 
0.080 ±0.011 ±0.003 
0.077 ± 0.013 ± 0.006 


Pt 


50 < 6» < 65 and 3 < z < 59 


65 < 6» < 80 and 6 < 2 < 59 


0.100 - 0.125 
0.125 - 0.150 
0.150 - 0.175 
0.175 - 0.200 
0.200 - 0.225 
0.225 - 0.250 
0.250 - 0.275 
0.275 - 0.300 


0.235 ± 0.028 ± 0.008 
0.227 ± 0.020 ± 0.017 
0.198 ±0.017 ±0.008 
0.145 ± 0.013 ± 0.008 
0.111 ±0.011 ±0.004 
0.109 ± 0.012 ± 0.004 
0.098 ± 0.013 ± 0.004 
0.100 ± 0.016 ± 0.015 


0.227 ± 0.029 ± 0.023 
0.214 ± 0.021 ± 0.007 
0.174 ±0.016 ±0.006 
0.157 ± 0.015 ± 0.010 
0.112 ±0.013 ±0.006 
0.094 ± 0.012 ± 0.006 
0.091 ± 0.015 ± 0.003 
0.137 ± 0.023 ±0.004 


PT 


80 < 9 < 95andll < 2 < 59 


95 < < 125 and 31 < z < 59 


0.100 - 0.125 
0.125 - 0.150 
0.150 - 0.175 
0.175 - 0.200 
0.200 - 0.225 
0.225 - 0.250 
0.250 - 0.275 
0.275 - 0.300 


0.326 ± 0.043 ±0.021 
0.214 ± 0.022 ±0.011 
0.168 ±0.017 ±0.008 
0.127 ± 0.015 ± 0.007 
0.112 ±0.014 ±0.004 
0.114 ±0.018 ±0.004 
0.114 ± 0.022 ±0.004 
0.153 ± 0.035 ± 0.008 


0.240 ± 0.027 ± 0.009 
0.201 ±0.020 ±0.011 
0.141 ± 0.018 ± 0.009 
0.129 ±0.018 ±0.006 
0.168 ±0.030 ±0.006 
0.099 ± 0.027 ± 0.023 
0.151 ±0.047 ±0.009 
0.227 ± 0.079 ± 0.021 
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5 PlON PRODUCTION BY NEUTRONS 

As shown schematically in Fig. [T| neutrons are part of the hadron cascade initiated by the 
interaction of the beam proton with a nucleus, in much the same way as protons. These neutrons 
are vitally important for the calculation of the tt^ content of the hadron cascade. While in 
the interactions of ~1 GeV/c protons with nuclei 7r+ secondaries strongly dominate over n~ 
secondaries, the opposite is the case for the interactions of ~1 GeV/c neutrons with nuclei. 

There are two issues relevant for the conventional flux from vr" decays. The first is the 
level of neutron production within the cascade, the second is the level of tt" production by 
neutrons. On both issues exist quite old yet relevant experimental data. 

As will be discussed in detail in Section [6l we base our estimations of the neutrino fluxes 
originating from the LSND 'beam stop' on two widely used Monte Carlo codes, Geant4 [18| 
and FLUKA [19]. This raises the question whether these codes reproduce correctly neutron 
production by protons on nuclei, and tt production by neutrons on nuclei. 

As for the first issue, neutron production by protons, we checked Geant4 and FLUKA pre- 
dictions against the results of an experiment pO| that measured neutron production by 1 GeV 
protons on 18 nuclei at polar angles of 4°, 7.5° and 11.3°. While we found agreement with the 
predictions of FLUKA within 15% — what we consider satisfactory — , we found that Geant4 
strongly underestimates neutron production by protons. Therefore, in certain simulation op- 
tions discussed in Section [6} we chose to replace Geant4 predictions of neutron production by 
protons by the respective FLUKA predictions. 

On the second issue, vr" production by neutrons, the results of a pertinent experiment [21 1 
are quantitatively in conflict with what Geant4 and FLUKA predict, although qualitatively both 
codes comprise pion production by neutrons. The situation is shown for the case of the copper 
nucleus in Fig. [v] We chose to ignore the small discrepancy in 7r+ production since the rela- 
tive contribution of neutrons to 7r+ production is small, and to subsume resulting differences 
between FLUKA and Geant4 predictions under our final systematic error estimate. However, 
the discrepancies in tt^ production between experiment and the Geant4 and FLUKA predic- 
tions are much too large to be ignored, given the much larger relative contribution of neutrons 
to TT production. Therefore, in certain simulation options discussed in Section |6} we chose 
to increase by a factor of two both the FLUKA and Geant4 predictions of tt production by 
neutrons. 

To our knowledge, LSND neglected tt production by neutrons which contributed to their 
underestimation of the conventional Vc flux originating from the 'beam stop'. 



6 Simulation of the neutrino fluxes in the LSND experiment 

An estimation of the neutrino fluxes from pion and muon decays in the LSND 'beam stop' is a 
major undertaking. It requires a detailed simulation of the geometry and material composition 
of the 'beam stop', of the characteristics of the hadron cascade after the interaction of incoming 
beam protons, and of all relevant physics processes in the interactions of particles with matter 
and in particle decays. 

The geometry and material composition of the LSND 'beam stop' as it was configured in 
the years 1993-1995, is depicted in Fig. 2 of Ref. [1] which is reproduced here in Fig. [8] 

The three main materials of the LSND 'beam stop' are the water of the primary target, 
the copper of the beam catcher, and the iron of the water target vessel and the shielding. Our 
simulation refers to the years 1993-1995 because afterwards the water target was replaced by 
a "close-packed high-Z" target (Section II.A in Ref. [|T|) which is not defined in the necessary 
detail. One major uncertainty in the simulation are the "Isotope Production Targets" that varied 



15 



• jt"^ (exp.) 
o jt" (exp.) 
— FLUKA 
■ - - Geant4 



50 100 150 

Polar angle [deg] 




Fig. 7: Comparison of inclusive cross-sections of 7r+ and vr production by 600 MeV neutrons on copper 
nuclei between experiment and the predictions by FLUKA and Geant4. 
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Isotope Production 
Targets 



Proton 
Beam 




Upper Perimeter of Containment i — i 
Box Shown in Phantom I — I '^2' 



PLAN VIEW, NEUTRINO SOURCE □ High Z Material 

Fig. 8: "The layout of the A6 beam stop, as it was configured for the 1993-1995 data taking." (Figure 
and its caption copied from Fig. 2 in Ref. |1|.) 



considerably during the data taking. We tried to simulate on the average the effect of these 
target plates that consisted in general of rather heavy materials, by three molybdenum plates 
with 22.2 mm thickness each, spaced by 22 cm. 

An intriguing feature of the LSND 'beam stop' is the relatively large open space downstream 
of the water target which permit more pions to decay and thus create a larger conventional z/g 
background than necessary. 

With a view to keeping the likelihood of programming mistakes as low as possible, we de- 
veloped two nearly independent Monte Carlo programs. One was based on the Geant4 Monte 
Carlo tool kit [TSl and focussed on a precise description of the geometry and material com- 
position of the LSND 'beam stop' to the best of our knowledge, but had limited versatility. 
The other was a specifically written stand-alone program, 'LSNDsim', that employed a simpli- 
fied 'beam stop' geometry with reduced granularity, was fast and flexible, and hence permitted 
high-statistics studies of different experimental configurations. 

The only connection between our two Monte Carlo programs were the double-differential 
inclusive cross-sections of secondary protons, neutrons and pions, by incoming protons, neu- 
trons and pions. The default values of these cross-sections were determined using the Geant4 
and FLUKA programs, and implemented as options in 'LSNDsim'. Within 'LSNDsim', the 
default values could be replaced by different ones. 

The simulations of all relevant physics processes in the interactions of particles with matter 
and in particle decays were entirely independent. 

The two programs agreed in their predictions of the conventional flux to better than 10%, 
for a closely similar experimental configuration. We concluded that the 'LSNDsim' results 
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can be trusted. All results presented below were obtained by running 'LSNDsim' with high 
statistics. 

'LSNDsim' provided results from the following six 'core' options: 

1. G4 option: Geant4 cross-sections were used; 

2. G4CDD option: Geant4 cross-sections were replaced by measured HARP-CDP cross- 
sections, Geant4 cross-sections of neutron production by protons were replaced by FLUKA 
cross-sections, and tt^ production by neutrons was enhanced by a factor of two; 

3. FL option: FLUKA cross-sections were used; 

4. FLCDD option: FLUKA cross-sections were replaced by measured HARP-CDP cross- 
sections, and vr^ production by neutrons was enhanced by a factor of two; 

5. LS option: cross-sections according to the LSND parametrization were used, with pions 
produced only in the primary interactions of beam protons; and 

6. FLLSO option: FLUKA cross-sections and pions from four generations were used, how- 
ever with neutrons switched off and with pions produced in the primary interactions of 
beam protons according to the LSND parametrization. 

Geant4 and FLUKA represent today's state-of-the-art in the simulation of hadronic pro- 
duction and hadronic interactions in matter. Differences between the results from these two 
programs permits insight into systematic uncertainties. 

The partial replacement of Geant4 and FLUKA cross-sections by measured values increases 
the reliability of predictions of the conventional Uc flux. The main problem in this endeavour is 
that measured data exist only in restricted kinematical domains, for specific beam momenta and 
for a limited number of target nuclei. Hence interpolations and extrapolations still could not be 
avoided. 

HARP-CDP pion production cross-sections from the interactions of 1.5 GeV protons could 
be directly used for the primary interactions of beam protons in water and copper. In those 
kinematic domains where data were available, the Geant4 and FLUKA cross-sections were re- 
placed. An average correction factor from the measured kinematical domain was applied to 
correct on the average the unmeasured kinematical domains. For the re-interactions of sec- 
ondary protons of lower momentum, the same dependence on momentum was used throughout 
as implemented in the Geant4 and FLUKA cross-sections, respectively. The production of vr^ 
by neutrons was enhanced by multiplying the respective Geant4 and FLUKA cross-sections by 
a factor of two. Pion production by tt"*" was not modified since the pertinent HARP-CDP data 
are taken with incoming 1.5 GeV/c 7r+ which is too high in momentum to help improving the 
relevant cross-sections at a few 100 MeV/c momentum. 



The measured 7r^/7r+ ratio from a 60 cm long water target reported in Section 4.2 was used 
in a fit that modified the inclusive cross-sections of 7r+ and vr^ production by 1.5 GeV/c protons 
in water, reported in Section|4| such as to reproduce closely the measured /n^ ratio. 

Cross-sections for the materials iron and molybdenum, not measured by HARP-CDP, were 
obtained by interpolation of respective measurements for copper, tantalum and lead. For the 
interpolation an A^^^ law was used as found appropriate for pion production by low-momentum 
protons (see, e.g., the discussion in Section 6 of Ref. [ITOl). 



Table 10 lists the materials of the 'beam stop' and the probability of the occurrence of the 



primary inelastic interaction of the beam proton. 



Table 1 1 gives a breakdown of what happens to 7r+ and vr in four generations after the 



primary inelastic interaction of the beam proton, for the FLCDD option of 'LSNDsim'. 
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Table 10: Materials of the 'beam stop' and the respective probabiUty of the primary inelastic interaction 
of the beam proton therein. 



Material 


Probability 


Copper (beam catcher) 


45% 


Water (target) 


27% 


Molybdenum (representing isotope production targets) 


22% 


Iron (representing the vessel and inserts of the water target) 


5% 


Aluminum (structural material) 


< 1% 


Air (free space) 


< 1% 



Table 1 1 : Breakdown of the relative abundances of vr^ and vr^ in four generations after the primary 
inelastic interaction of the beam proton, and of what happens to them, for the FLCDD option of 'LSND- 
sim'. 







71+ 


TT 


1st generation 




77% of total 


47% of total 




Inelastic interaction 




41% 


34% 




Decay in flight 




2% 


2% 




Decay at rest (tt"*") or disappear (n' 




57% 


64% 


2nd generation 




20% of total 


40% of total 




Inelastic interaction 




32% 


28% 




Decay in flight 




2% 


1% 




Decay at rest (tt^) or disappear (tt" 


') 


66% 


71% 


3rd generation 




3% of total 


11% of total 




Inelastic interaction 




25% 


26% 




Decay in flight 




2% 


1% 




Decay at rest (tt"*") or disappear (tt" 


-) 


73% 


74% 


4th generation 




<1% of total 


2% of total 




Inelastic interaction 




22% 


24% 




Decay in flight 




1% 


1% 




Decay at rest (7r+) or disappear (tt" 


-) 


77% 


75% 
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Negative muons that come to rest are either captured or decay in orbit. The capture prob- 
ability strongly depends on the atomic number Z but is generally well known. We took the 
capture probabilities Acapt from Ref. [22|. The decay-at-rest probabilities of negative muons 
Adec/(Adec + Acapt) that wc uscd iu our simulation, are listed in Table 12 in decreasing order. 
For the estimation of the conventional flux, it is important that the relative weight of each 



Table 12: Decay-at-rest probabilities Adec/(Adoc + Acapt) of negative muons in materials of the LSND 
'beam stop', as used in the 'LSNDsim' program. 



Material 


Decay-at-rest probability 


Air 


0.862 


Water 


0.818 


Aluminum 


0.396 


Iron 


0.094 


Copper 


0.074 


Molybdenum 


0.047 



material in which negative muons come to rest, is correctly simulated. The effect on different 
generations will be different, for their different spatial extensions and hence different sampling 
of materials. Table 13 gives the fractions of the conventional flux that originate from /x^ 
decays in different materials, and the fractions that originate from the decays of /i^ of different 
generations, for the FLCDD option of 'LSNDsim'. 



Table 13: Fractions of the conventional flux that originate from decays in different materials, and 
in different generations, for the FLCDD option of 'LSNDsim'. 





Fraction of conventional 


Material 






L-on 


51% 




Copper 


29% 




Water 


12% 




Molybdenum 


4% 




Aluminum 


4% 




Air 


<1% 


Generation 






1st 


74% 




2nd 


22% 




3rd 


4% 




4th 


<1% 



In the following, we shall undertake the instructive comparison of the predictions for 7r+ and 
TT production in the four options G4CDD, FLCDD, LS and FLLSO of 'LSNDsim'. We note 
that despite imposing HARP-CDP cross-section results in the G4CDD and FLCDD options, 
there are enough differences in the hadron production characteristics in Geant4 and FLUKA 
that one cannot expect the same results from G4CDD and FLCDD. These are differences in the 
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double-differential cross-sections of hadron production, especially by vr^ and neutrons where 
HARP-CDP have no data to contribute, and their momentum dependence; and differences in 
cross-sections in kinematic regions that are outside of the regions accessible by the HARP 
experiment. 

Figure |9] demonstrates the contributions from different generations to the overall pion pro- 
duction in the LSND 'beam stop', for the example of the FLCDD option. We recall that 'LSND- 
sim' employs four generations of hadron production and Fig. |9] makes clear, especially for tt , 
why this was deemed necessary. 

Figure 10 addresses the issue of differences in the overall yields (summed over four genera- 
tions) of pion production in the LSND 'beam stop' between the four options G4CDD, FLCDD, 
LS and FLLSO. 

While the differences in overall 7r+ production are at the 20% level, the differences in over- 
all TT production between G4CDD and FLCDD on the one hand, and LS and FLLSO on the 
other hand, amounts to a factor of two. The cause is mainly the neglect of pion production by 
secondary neutrons in the LS and FLLSO options. At first sight this difference looks dramatic 
for the conventional flux. However, it has less dramatic consequences when taking into ac- 
count that pions that are produced by neutrons cannot occur in the first two generations, and 
hence such pions miss the chance of decaying in the ample free space after the water target 
(c.f. Fig. [8]). Instead, their decay path is shortened to the interaction length in the copper beam 
catcher, which reduces the impact of tt" produced by neutrons since arise solely from 7r~ 
that decay in flight before coming to rest. 

Table [14] presents the salient results of our simulation work and a comparison with the 
results from the 1993-1995 data taking of the LSND experiment, as stated in Ref. yj- It lists 
the absolute rates of 7r+ and vr^, the ratio 7r^/7r+, the absolute rates of z/^j from vr"*", and of Ufj, 
from TT^, that decay in flight. The z/^ and cases hold for an energy above 123.7 MeV and 
113.1 MeV that are thresholds for the reactions z/^ -i- ^^C fi^ + ^^N* and z/^ -i- p — )■ /i+ -i- 
n, respectively, that lead to backgrounds as will be discussed in Section |7] Neutrinos of such 
energy can only arise from pions that decay in flight (termed 'DIF' by LSND). Further the 
absolute rates, from zero to the maximum energy of 52.8 MeV, of z/^ and z/c, and of z/^ and 
z>e, are given that originate from fi^ and yU^, respectively, that decay at rest (termed 'DAR' by 
LSND). We recall that the rate from /i+ [DAR] decays is rather safe to calculate as it depends 
essentially on the integral of all 7r+ produced in the hadronic cascade in the 'beam stop'. By 
contrast, the /i^ [DAR] rate (that is the centrally important one), is rather uncertain to calculate, 
as it requires a vr^ decaying in flight into /x" -i- z/^, and this /i^ decaying at rest into e~ + Uc + 
in competition with much stronger nuclear capture. The ratio z/o/z?/^ from muon [DAR] is given 
to emphasize the smallness of the conventional z/c flux. 

As for the rationale of the comparison with LSND results we stress the following. We 
cannot claim that our simulation program describes the LSND 'beam stop' in sufficient detail, 
for we rely solely on information that has been published. Therefore, a comparison of absolutely 
normalized rates from our simulation with absolutely normalized rates published by LSND is of 
limited interest. Also, in the centre of interest is not the detailed geometry of the LSND 'beam 
stop' and the detailed treatment of particle transport, but rather the hadronic cascade initiated 
by the primary interaction of the beam proton. Therefore, we replace the absolute comparison 
by a relative comparison where details of the 'beam stop' geometry and of particle transport 
cancel in first approximation. 

We claim that the 'LSNDsim' options LS and FLLSO encompass the hadronic physics model 
that was employed by LSND. We also claim that the G4CDD and FLCDD options of 'LSND- 
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Fig. 9: Absolutely normalized yields of vr"*" (upper panel) and vr (lower panel) production of the first 
four generations after the initial interaction of the beam proton, for the FLCDD option. 
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Fig. 10: Absolutely normalized yields of overall 7r+ (upper panel) and tt (lower panel) production as 
predicted by the four options G4CDD,FLCDD, LS and FLLSO of 'LSNDsim'. 
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sim' represent more closely the hadronic cascade that developed in the LSND 'beam stop'. 
Our argumentation involves four steps: 

1. we demonstrate that the LS and FLSSO options reproduce reasonably well the results 
published by LSND (it is not required that they reproduce exactly the LSND results); 

2. we take the average of LS and FLSSO results as reference 'LSND results'; 

3. we take the average of G4CDD and FLCDD results as the best estimate of what LSND 
should have measured; 

4. we apply the difference between the average of G4CDD and FLCDD, and the average 
of LS and FLLSO, to the published results of LSND and consider this as the result that 
LSND should have obtained. 

In Table [15| we compare the absolutely normalized rates of pions and neutrinos from the 
1993-1995 data taking of the LSND experiment as published in Ref. [[T|, with the averages of 
the 'LSNDsim' options LS and FLLSO ('Emulation'), and G4CDD and FLCDD ('Best esti- 
mate'). 

We note the increase by a factor of 1.6 of the conventional rate of the 'best estimate', 
0.919 X 10"^^ (pot • cm^)-!, over the one of the 'emulation', 0.585 x 10"^^ (pot ■ cm^y^. 

In our simulation, the ratio 7r^/7r+ refers to the average pion production in the LSND 'beam 
stop'. Our ratios in the LS and FLLSO options are around 0.2, in conflict with the ratio 0.12 
that is repeatedly quoted by LSND (e.g., in Table VII of Ref. [ 1 1). We conjecture that the latter 
ratio must refer to something different than the average pion production in the entire 'beam 



stop' — that is why we put in brackets the respective LSND number in Table 15 



Table 16 lists the uncertainty estimates that we felt appropriate to assign to our 'LSNDsim' 
simulation. 

To the systematic errors listed in Table 16, a 5% error for the 7r+ chain and a 9% error for 
the TT chain from HARP-CDP data statistics is to be added. 

Adding all errors in quadrature, the overall error for neutrinos from the 7r+ chain is 15%, 
and 24% for the vr" chain. These are the errors that we consider appropriate for our prediction 
of absolute integral neutrino fluxes from the LSND 'beam stop'. 

Not all listed errors are relevant when one refrains from absolute predictions and asks what 
error to be assigned to the relative comparison between 'LSNDsim' options LS and FLLSO on 
the one hand, and the LSND Monte Carlo on the other hand. We consider that the errors are for 
the 7r+ chain not less than 9%, and for the tt^ chain not less than 11%, respectively. 
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Table 14: Absolutely normalized rates of pions and neutrinos from the 1993-1995 data taking of the 
LSND experiment from several 'LSNDsim' options; the rates are given per 'pot' (proton-on-target) and 
per unit of area (cm^); DDF rates refer to pion decay in flight, with the cut E^, > 123.7 (113.1) MeV 
applied; DAR rates refer to muon decay at rest, with no cut in E,y appUed. 





LS 

("Emulation A") 


FLLSO 
("Emulation B") 


7r+ Inot^^l 

TT" [pOt"^ 


0.096 
0.020 


0.121 
0.024 




0.206 


0.195 


7r+ [DIF] z/„ [10-^^ (pot-cm^)-!] 
71- [DIF] z/^ [10-12 (pot-cm^)-!] 


2.78 
2.28 


3.10 
3.08 


/x+ [DAR] z7^, z/e [10-9 (pot-cm2)-i] 
a- [DAR] Pe [10-^2 (pot-cm^)-^] 


0.520 
0.547 

10.5 X 10-^ 


0.673 
0.622 

9.24 X 10-4 




G4 


G4CDD 
("Best estimate") 


7r+ lnot~^l 

TT" [por^] 


0.152 
0.030 


0.140 
0.050 




0.196 


0.358 


TT"^ [DIF] z/^ [10"^^ (pot-cm^)-!] 
TT" [DIF] z/^ [10-12 (pot-cm^)-!] 


6.77 
5.61 


5.57 
7.83 


/i+ [DAR] z7^, z/e [10-9 (pot-cm2)-i] 
11" [DAR] z/^, Pe [10-12 (pot-cm2)-i] 


0.834 
0.675 

8.09 X 10-4 


0.777 
0.956 
12.3 X 10-4 




FL 


FLCDD 
("Best estimate") 


TT"^ [pOt-l] 
TT" [pOt-l] 


0.141 
0.030 


0.130 
0.045 


77-/77"*" 


0.213 


0.343 


77+ [DIF] z/^ [10-11 (pot-cm2)-i] 
77- [DIF] P^ [10-12 (pot-cm2)-i] 


2.08 
2.23 


1.83 
2.44 


//"^ [DAR] P^, i/e [10-9 (pot-cm2)-i] 
^x- [DAR] i/^, Pe [10-12 (pot-cm2)-i] 


0.833 
0.642 
7.71 X 10-4 


0.763 
0.881 
11.5 X 10-4 
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Table 15: Absolutely normalized rates of pions and neutrinos from the 1993-1995 data taking of the 
LSND experiment as published in Ref. [1], and from the average of the 'LSNDsim' options LS and 
FLLSO ("LSND emulation"), and G4CDD and FLCDD ("Best estimate"); DIF rates refer to pion decay 
in flight, with the cut > 123.7 (113.1) MeV applied; DAR rates refer to muon decay at rest, with no 
cut in applied. 





LSND 
published 
(1993-1995) 


LS / FLLSO 
("Emulation") 


G4CDD / FLCDD 
("Best estimate") 


7r+ [por^] 
TT" [por^j 




0.109 
0.022 


0.135 
0.048 




(0.12) 


0.203 


0.352 


n+ [DIF] [10^11 (pot-cm^)-!] 
TT- [DIF] [10-12 (pot-cm^)-!] 


1.48 
1.57 


2.94 
2.68 


3.70 
5.14 


/i+ [DAR] z7^, [10-9 (pot-cm2)-i] 
/i- [DAR] z/^, [10-12 (pot-cm2)-i] 


0.8 
0.65 

8.1 X 10-^ 


0.597 
0.585 
9.80 X 10-^ 


0.770 
0.919 

11.9 X 10-^ 



Table 16: Estimated uncertainties in the 'LSNDsim' simulation. 



For both the 7r+ and tt chains 


HARP-CDP systematics 


8% 


Rescaling of unmeasured regions 


3% 


'Beam stop' geometry 


3% 


Cross-section binning 


6% 


Cross- section interpolation 


6% 


Hadron multiplicity 


5% 


Extrapolations to different nuclei 


4% 


No 5th hadron generation 


1% 


In addition for the tt chain 


'LSNDsim' statistics 


2% 


TT^ production by neutrons 


3% 


neutron production by protons 


3% 


muon capture 


4% 


Decay space geometry 


5% 


Momentum spectrum 


15% 



26 



7 New estimates of Backgrounds 

The 'LSND anomaly' is the synonym of a rate of z/c + P ^ + n events that is much larger 
than expected. In such a situation, a thorough discussion of backgrounds is mandatory. 

The reaction Uc + p ^ + ^ has two types of backgrounds. Background I is caused by 
the flux of conventional origin. According to the discussion on the conventional z/g rate in 
Section [6[ this background is to be increased by a factor of 1.6 and its relative error increases 
from 20% to 29%. 

Background II consists primarily of events with low-momentum muons in the final state that 
are 'invisible' (we remind of the copious fluxes of and z/^ in the LSND experiment that make 
this background important; normally, a muon is identified as such in the LSND detector and the 
respective event does not enter the sample of signal candidates). LSND state several reasons for 
a muon to be 'invisible'. The by far most important one is a kinetic muon energy < 3 MeV 
(Section VILA in Ref. [ 1 ]). Of prime interest is the reaction z/^ + p — > /i+ + n since it leads to 
a neutron in the final state. Less important are the reactions + ^^C — )■ /i"^ + ^^B* and z/^ + 
^^C + ^^N*. The common feature is that the final-state muon is 'invisible', that final-state 

particles (or decay products from final-state nuclei) lead to or fake the neutron capture process, 
and that the event is triggered by an electron from muon decay. 

The difference between LSND's and our assessment of Background II event numbers is less 
in the expectation values but rather in their uncertainties. 

Only z/^ with energy above 123.7 MeV can initiate the reaction z/^ -i- ^^C + ^^N*. The 

thresholds for the reactions z/„ -i- p — )• /i+ -i- n and z/„ -i- ^^C + ^^B* are 113.1 MeV and 



1 19.7 MeV, respectively [23 1. That means that only the high-momentum end of the 7r+ and it~ 
spectra can give rise to such neutrinos. These small portions of the 7r+ and vr^ spectra have a 
larger uncertainty than the integral of these spectra, though. 



This is illustrated in Fig. 1 1 which shows for the FLCDD option of 'LSNDsim' on a log- 
arithmic scale the whole tt spectrum and the portion of the spectrum that is cut out by the 
requirement of an accepted z/^ from tt" [DIP] with energy greater than 1 13.1 MeV^^ The large 



uncertainty of this spectrum is highlighted in Fig. 12 which shows the differences between the 
G4CDD, FLCDD, LS and FLLSO options of 'LSNDsim' . The uncertainty is a priori at the 
100% level^^ Although we apply eventually an error smaller than 100%, for compliance with 
LSND's 'constraints' on neutrino fluxes that will be discussed in Section[8} we consider LSND's 
claims unconvincing that the numbers of 7r+ and tt^ that decay in flight and give rise to muons 
with < 3 MeV are known with an error as small as 15%. 

LSND calculate Background II as a small fraction of all events above the respective neutrino 
energy thresholds. We consider that this approach renders Background II more dependent on 
the detailed shape of the high-energy ends of the pion spectra than necessary. We consider it 
safer to calculate Background II from the small portion of the neutrino spectrum between the 
neutrino energy threshold and 4 MeV above^^ . 



^^We show the tt^ spectrum because the contribution of the reaction i/^ + p — > /i+ + n to Background II is 
particularly large. 

^^The argument that one cannot conclude from the spectrum integral on the size of the small high-energy portion 
of the pion spectrum can also be turned around: one cannot conclude from the high-energy portion on the size of 
the spectrum integral. 

^^The value of 4 MeV is approximate and beset with instrumental uncertainties; we bypass these uncertainties by 
considering not absolute numbers but relative numbers; we determine a factor by dividing the average of G4CDD 
and FLCDD predictions by the average of LS and FLLSO predictions, and apply this factor as a correction to 
LSND's estimate of Background II. 
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Fig. 11: The 'LSNDsim' prediction (FLCDD option) of the whole tt spectrum and of the portion of the 
spectrum that is cut out by the requirement of an accepted from tt" [DIF] with energy greater than 
113.1 MeV. 
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Fig. 12: Predictions of the energy spectra of accepted i/^ from tt [DIF] with energy greater than 
113.1 MeV, in the G4CDD, FLCDD, LS and FLLSO options of 'LSNDsim'. 
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In Table 17 we summarize the reactions and respective event numbers that constitute Back- 
ground II, and compare LSND's claims with our assessments. 



Table 17: LSND Background II event numbers 



Reaction 


Background II type 


No. events 




LSND publ. 


This paper 


P ^^/^^ n 


< 3 MeV 


8.2 


10.8 ±8.0 




i^N*, < 3 MeV, 


1.4 


1.8 ±1.8 




capture 




0.2 ±0.2 




< 3 MeV 


0.4 


0.5 ±0.5 


Otherwise missed muon 




0.4 ±0.14 


0.4 ±0.14 


/i~ e" h>c 1^,1, vr" e" 


z/e events 


0.1 ±0.1 


0.1 ±0.1 


Sum 




10.5 ±4.6 


13.8 ±8.2 



For completeness. Table [TT] lists in addition to reactions with misidentified low-momentum 
muons two more — minor and uncontroversial — contributions to Background II. 

In Table 18 we compare Backgrounds I and II as published by LSND in Ref. [HI, with the 
respective results obtained from 'LSNDsim'. The 'LSND anomaly' which is LSND's measured 
'beam excess' of 117.9±22.4 events after the subtraction of Backgrounds I and II, changes from 
the published value of 87.9 ± 23.2 to our value of 73.5 ± 25.4. Its significance reduces from 
3.8 a to 2.9 a. 



Table 18: Backgrounds I and II as published by LSND and as obtained from our 'LSNDsim' program. 





LSND published 


This paper 


'Beam excess' 


117.9 ±22.4 


Background I 


19.5 ±3.9 


30.6 ±8.8 


Background II 


10.5 ±4.6 


13.8 ±8.2 


'LSND anomaly' 
Significance 


87.9 ±23.2 
3.8 cT 


73.5 ±25.4 
2.9 cr 



We reiterate that in a subsequent paper Q we shall argue that the 'LSND anomaly' and its 
significance is even smaller. 

8 On LSND's 'constraints' on neutrino fluxes 

LSND claim that their neutrino fluxes are correct within the errors quoted in their final physics 
publication [1| because certain experimental 'constraints' are met. They list these 'constraints' 



in Table IV in Ref. [1] that is reproduced here as Table 19 However, there is no quantitative 
follow-up of the 'constraints'. 

In the following, we shall argue that (i) our estimate of Background I from the interactions of 
conventional u^, although larger by a factor of 1.6, is not in conflict with LSND's 'constraints'; 
and that (ii) LSND's 'constraints' concern primarily not Background I but the less important 
Background II from misidentified z/^ and z/^ interactions. 
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Table 19: "Cross section uncertainties for the neutrino reactions with two-body final states that occur in 
LSND. The cross sections for these processes are known accurately because either related measurements 
can be used to constrain the matrix elements or only fundamental particles are observed. Also shown are 
the corresponding neutrino flux constraints." (Table caption and table contents copied from Table IV in 
Ref. Q). 



Process 


a Constraint 


a Uncertainty 


Flux Constraint 




Standard model process 


1% 


/i+ z/cZ/^e+ DAR 


i^C (z/e,e-)i2Ng, 




5% 


/i+ — )■ z/eZ/^e+ DAR 




12n 


5% 


71+ -> z/^/i+ DIF 


P (z^M'/^") n 


neutron decay 


5% 


7i~ — 7- I'fifi" DIF 



These are the cornerstones of our argumentation: 

1. On the Vc ^^^d i?^ fluxes from fi+ [DAR] decays 

LSND claim in Ref. [24j that from the observed number of events from the exclusive 
reaction z/g + ^^C ^ e~ + ^^Ngs, the total n+ flux is known to 11% through the chain 7r+ 

We are not disputing this 11% claim. However, we note that LSND claim with 7% a 
better precision than 11% on the neutrino flux from fi+ [DAR] decays, and therefore on 
the total z/c and z/^ fluxes. The better precision does not stem from the said 'constraint' 
but from the calibration experiment E866 [25| (which we dispute, see Section |9]). 
LSND claim in Table IV in Ref. |[T| also the elastic scattering process z/ + e — )■ z/ + e~ 
to constrain the Ue and fluxes from /i"*" [DAR] decays. We hold that the statistics are 
scarce and translate into a constraint on the total 7r+ flux at the 20% level — which brings 
nothing new. 

2. On the z/^ and u^i fluxes from n+ [DIF] and it~[DIF] decays, respectively 

LSND claim in Ref. [26 1 56.8 ± 9.6 (17%) events from the exclusive reaction z/^ + ^^C 
H~ + ^^Ngs and a measured cross-section of (6.6 ± 1.0 ± 1.0) x 10""^^ cm^ (21%); for the 
determination of this cross-section, the z/^ flux above the muon threshold of 123.7 MeV 
is assumed by LSND to be known to 15%. 

We hold that all these z/^ stem from n+ [DIF] which is only 2% of all n+ and the increase 
from the claimed 7% precision on the number of all 7r+ to the claimed 15% precision on 
the number of z/^ at the high-energy end above 123.7 MeV from n+ [DIF] decays is too 
optimistic. We dispute LSND's claim to have measured the cross-section of z/^ -i- ^^C — )■ 
^- + i^Ng^ with 21% precision. 

We turn the argument around and accept that the measured number of events from the 
exclusive reaction z/^ -i- ^^C fi^ + ^^Ngs together with its theoretically well-known 
cross-section (5%) constrains LSND's z/^ flux above 123.7 MeV to 17% precision. How- 
ever, we note that this precision refers to all z/^ above 123.7 MeV but does not constrain 
in the same way the z/^ spectral shape as a function of energy; specifically, the small band 
between the muon threshold and 4 MeV above that is relevant for the estimate of Back- 
ground II must have a much larger uncertainty than 17%. Before we discuss how much 
larger, we consider the 'constraint' on the flux. 

LSND measured (see Table V in Ref. p6| ) in the 1993-1995 exposure the sum of all 
events from the reaction + p ^ ii+ + n and the two inclusive reactions z/^ -i- ^^C — )■ n~ 
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X and + ^^C — )■ yU+ X as 1924 ± 47 events, where — according to cross- sections from 
theory since no experimental distinction is possible — 140 events come from the reaction 
z/^ + p — J- /i+ + n, 46 events from the reaction z/^ + ^^C — )• /i+ + X and 1738 events from 
the reaction +^^C n~ + X. With a 17% precision on the flux above the muon 
threshold of 123.7 MeV and with an estimated fraction of all events of 10.7% due to z/^ 
reactions, no constraint can be derived on the z/^ flux. 

Yet LSND claim another constraint from the measured number of events with a muon 
together with a correlated neutron from the three reactions z/^ + p — )■ yU+ + n, z/^ + ^^C — t- 
yU+ + n + X and z/^+ ^^C — )• + n + X. This approach is motivated by the (theoretical) 
expectation that the requirement of a correlated neutron much reduces the preponderance 
of the z/^ reaction with ^^C. The observed number of events is 210 ± 35 events in the 



1993-1995 exposure (see Table VI in Ref. |26|). The calculated repartition of events 
among the three reactions is 210 = 140 + 36 + 34. LSND conclude from this that "The 
observed number of events with neutrons also rules out a z/^ flux much bigger than that 
calculated by the beam Monte Carlo simulation" (Section VII in Ref. [j26|). 
We observe first that LSND do not state numerically what 'much bigger' means; second, 
the argument depends on whether the fraction of events with a correlated neutron is de- 
termined correctly (which we dispute); and third, the argument hinges on the correctness 
of the assumed fraction of correlated neutrons (which is certain as 100% for z/^ -i- p — )■ /x"*" 
-I- n, while predicted by theory as 60% for z/^ -i- ^^C — )■ /i+ -i- X and as 6% for z/^ -i- ^^C — 
fi-+X). 

With an uncertainty of the z/^ flux above 123.7 MeV of 17%, with an uncertainty of the 
z/^ -I- ^^C — )■ /i^ -I- X cross-section of 100%, and with neutron branching fractions of 
(60 ± 20)% and (6 ± 4)%, respectively, there follows an uncertainty of 32% of the z/^ flux 
above 113.1 MeV. 

The z/^ and z/^ fluxes are only relevant for the estimation of Background II from 'invisible' 
muons; because this background stems only from the narrow band of neutrino energies 
between threshold and 4 MeV above threshold, and the error on the total fluxes above 
threshold says little about the systematic errors of event numbers from neutrinos in the 
narrow bands, we double and then round the errors for the fluxes in the narrow bands: 
35% for z/^ and 60% for z/^. These errors are used to calculate the error of Background II 
listed in Tables [IT] and [H 

The most important conclusion is that there is no strong constraint on the flux from /x^ 
[DAR]. The 32% uncertainty of the z/^ flux above the 113.1 MeV threshold is well consistent 
with our claims of the increase by a factor of 1.6 of the conventional z^o rate (that stems from 
the integral of all tt and not from the high-momentum end of tt decaying in flight), and of 
the increase of its relative error from 20% to 29%. 

While LSND's 'constraints' have little impact on the all-important conventional rate, they 
concern primarily the rate of background events with an 'invisible' muon in the final state, and 
its error. 



9 Critique of LSND's neutrino flux calculations 

Our critique of LSND's neutrino flux calculations focusses on two aspects. The first is the 
inadequate handling of contributions from higher-generation particles in the hadronic cascade 
after the primary interaction of beam protons in the LSND 'beam stop' , primarily from neutrons. 
The second is the assignment of too optimistic systematic errors of the calculated neutrino fluxes 
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(We quote: "Calculations of /i+ DAR fluxes are uncertain at the 7% level, while vr^ DIF fluxes 
and ji^ DAR fluxes are uncertain to 15%", c.f. Section II.D in Ref. yj). 

The experimental pion production cross-sections that were used as input to LSND's Monte 
Carlo neutrino flux calculations [[T4l[T5| were thin-target cross-sections. However, the LSND 
'beam stop' represented a thick-target geometry. How was this important difference taken into 
account? 

We quote again: "The proton beam degradation within the target is described by the inelas- 
tic proton cross section, for which we use the energy-independent cross section measured in 
neutron-nucleus collisions... Each inelastic proton collision is assumed to reduce the energy of 
the proton by an amount appropriate for pion production... Protons are followed, losing energy 
either through colli sional or ionization energy loss, until reaching the minimum energy for pion 
production... Pion production through secondary proton interactions accounts for approximately 
10% of the total." (Section 2.2 in Ref. y4J). 

There was clearly awareness that secondary protons play a role and this was taken into ac- 
count in the way described. However, it seems that LSND underestimated the importance of 
hadrons up to the fourth generation, and the importance of neutrons especially for tt^ produc- 
tion. 

LSND claim that their Monte Carlo program was tuned to reproduce pertinent experimental 
data from the E866 experiment [25 1, and use this as justification of their error assignments. 

We hold that the E866 experiment was sensitive only to the integral of 7r+ production which 
is rather uncontroversial. E866 claimed a 7.6% precision on this integral. We note that the 
LSND Monte Carlo program, when used for a simulation of the E866 experiment, "was un- 
able to reduce the predicted leakage from the side without making unacceptable changes to the 
pion decay inside the stack" (Section 3 in Ref. |[14|). The discrepancy between observed and 
predicted pion leakage, apparent in Fig. 1 1 of Ref. [[T4| suggests that even after adjustment of 
overall normalization, the LSND Monte Carlo exhibits serious shortcomings in the predictions 
of the momentum dependence of pion spectra. As a consequence, the systematic precision of 
the flux from stopped 7r+ must be worse than the 7.6% precision claimed by E866. 

We hold further that the E866 experiment is hardly sensitive to the vr component of the 
hadronic cascade, so there is from E866 no experimental check of LSND's claims on the amount 
and the uncertainty of tt production in the LSND 'beam stop'. 



10 Conclusion 

For the importance of the 'LSND anomaly' that is in stark conflict with the Standard Model, we 
have undertaken a re-evaluation of the data analysis of the LSND experiment. We compare an 
emulation of this experiment where we try to reproduce closely the LSND data analysis, with 
results of a simulation where the model of hadron production is considerably more realistic. 
We use new experimental cross-sections of pion production by protons from HARP-CDP, and 
four decades old data on pion production by neutrons. The more realistic model of hadron 
production leads to a considerably larger conventional Ue flux than was estimated by LSND. 
Taking merely into account the more realistic model of hadron production, the significance of 
the 'LSND anomaly' reduces from 3.8 a to 2.9 a. 
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